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Abstract 

It is shown that recent results from the Super-Kamiokande detector 
constrain the tau neutrino diagonal magnetic moment to u Ut < 1.3 x 
10~ 7 ub for the case of — > v T interpretation of the atmospheric neutrino 
anomaly. It is pointed out that the large magnetic moment of the tau 
neutrino could affect further understanding of the origin of the anomaly. 



Recent results from the Super-Kamiokande (S-K) detector give evidence for 
neutrino oscillations |I] - @. One of the favourable interpretations of the atmo- 
spheric neutrino anomaly suggested by the S-K collaboration is related to the 
existence of — > v T neutrino oscillations. This implies that tau neutrinos have 
non-zero masses and therefore may also have a non-zero diagonal magnetic mo- 
ment ( see e.g. ref. ||). As a consequence, massive tau neutrinos could manifest 
themselves in terrestrial experiments through the effect of —>■ v T neutrino oscil- 
lations, and, if the magnetic moment value is large enough, through tau neutrino 
electromagnetic interactions. 

In ref. (see also [0), it was shown that the combined existence of — > v T 
(and/or v e — > v T ) oscillations and non-zero magnetic moment of the tau neutrino 
would increase the total rate of events in v^{v e ) neutrino- electron scattering 
experiments. It was used to constrain the mixing angles of the tau neutrino with 
neutrinos of another flavour. In this Letter it is shown that this effect can also 
be used to constrain the magnetic moment of the tau neutrino from the S-K 
atmospheric neutrino data. 

Assuming that a muon neutrino beam has a component of tau neutrinos due 
to — ► v T oscillations, in the case of two-neutrino mixing, neutrino states evolve 
with time t as 



\u{t) >= a(t)\^ > +6(t)|z/ r ,/v + > (1) 
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where > and \v T > denote weak eigenstates of and v T neutrinos, and 
a 2 (t), b 2 (t) are the probabilities of finding or v T in the beam at a given moment 
t, respectively. It is assumed that a 2 (0) = 1 at t — 0. The probability b 2 (t) depends 
on the parameters of — v T oscillations as [§]: 

b 2 (t ~ -) = P(z/ M -> v T ) = sinHd^sin 2 ^^- (2) 

or 

2 „ 2 1.27Am 2 (e1/ 2 )L(A;m) 
- i/ T ) « S m 2 2VW Eu (GeV) (3) 

where sin 2 26 flT is the mixing angle, Am 2 = |m 2 — m|| is the difference of the 
squares of the mass eigenstates in eV 2 , E v is the neutrino energy in GeV, and 
L is the mean distance between the neutrino source and the detector in km. In 
the above formula it is also assumed that magnetic field B is weak enough not 
to affect the probability of oscillations Piy^ —>■ f T ), i.e. |§: 



Am 2 /2E U > fi UT B (4) 

Then if the magnetic moment of the v T exists, it will contribute to a non- 
coherent part of the v T e~ scattering cross section via the reaction that changes 
the helicity of the tau neutrino (hence right-handed neutrino states should exist). 
This might result in a contribution to the observed deviations from unity of 
the flavour ratio of the atmospheric neutrino flux, R = ((N v + N v )/(N Ue + 
Nu e )) da J((N u ^ + N Vf )/{N Ve + N Ve )) MC , in the S-K detector for both data and 
Monte Carlo (MC) simulation. Here, N v jV and N Ve)Ve are the number of muon 
-like (//-like) and electron-like (e-like) fully-contained events in the S-K detector 
I- 

Indeed, since the electromagnetic cross section is orders of magnitude larger 
than the weak cross section, even a small fraction of tau neutrinos with non-zero 
magnetic moment in the atmospheric neutrino flux could lead to an observable 
excess of isolated electrons in atmospheric neutrino interactions in the S-K detec- 
tor whose signature is identical to that of e-like events. Note that the magnetic 
moments of muon or electron neutrino are experimentally proved to be too small 
to give a noticeable contribution to the neutrino interaction rate. The produc- 
tion rate of isolated electrons via v T e~ scattering in the detector depends on 
the probability Piy^ — > v T ) of finding v T neutrinos in the atmospheric neutrino 
flux. This probability can be calculated from the neutrino survival and transition 
probabilities using Eq.(2). 

The neutrino-electron scattering process via magnetic moment has a cross 
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where E e is the electron energy, fi u is the neutrino magnetic moment and 
fiB = e/2m e is Bohr magneton. For the case of constant E v the integral cross 
section is 
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are the high and low electron energy cuts, respectively. 



Note that the integral cross- section <r M depends very weakly on the neutrino 
energy. It rises only logarithmically with the neutrino energy, while the total 
neutrino cross section rises linearly with E v . Thus, it is advantageous to search for 
a tau neutrino magnetic moment using the low energy (sub-GeV, E visi bi e < 1-33 
GeV) S-K data 0. 

Approximately, all atmospheric neutrino giving rise to the fully contained [i- 
like events in the S-K detector have energy above 400 MeV, see e.g. Ref. [ 1 1 



(MC studies showed the mean neutrino energy for CC interactions to be about 
800 MeV for /i-like events |jlj|). Thus, for the lowest electron energies analysed in 
the S-K, from E™ in = 100 MeV to E™ ax = 200 MeV (first bin of the histogram 
in Fig. 4 (a) corresponding to e-like events, ref. JI]), the cross section of v T — e 
scattering due to non-zero magnetic moment is 



a M = 1.2 x 10- 39 (fi u /10- 7 fi B ) 2 cm 2 (7) 

Here, we assume the tau neutrino energy E v = 400 MeV, thus the integral 
cross section is underestimated. 

The total number AN e of e-like events from v T — e scattering in the first energy 
of Fig.4(a) in |IJ can be written in the form 

AN e = J k ■ f v ■ a M ■ e ■ dn ■ dE u = hN^a^ (8) 

where ko is a factor related to the ^-target mass, f u is tau neutrino flux, 
the cross section is constant, e is detection efficiency which is practically 
independent of energy for E e > 100 MeV k\ is a factor corresponding to 
the convolution of detector acceptance, detection efficiency and //-target mass 
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and N Vt is the total number of tau-neutrinos crossing fiducial volume of the S-K 
detector. 

The number of tau-neutrinos is estimated using the number of the initially 
produced muon atmospheric neutrino: 



N Ut = (N^+N v JP(u,^u t ) (9) 

The observation of the small value of R = 0.61 ± 0.03(stat.) ± 0.05(syst.) |l| 
and zenith angle dependent deficit of //-like events suggests, that in the case of 
two -neutrino oscillation scenario of z/ M — > v T more than 30% (90%C.L.) of the 
initially produced muon neutrinos arrive at the S-K detector as tau neutrinos 
IfJ. |3|. We assume that the average probability of oscillations P{y^ — > v T ) = 0.3. 

Muon-neutrino quasi-elastic events detected in the S-K detector were used to 
estimate the expected total number of muon neutrino. Similarly to Eq.(8) one 
can write 



NT = UN»y q . e . + NvtfJ (l - P(i/„ -> v T j) 
= k 2 (N u ^ + N Vfi ) (l - P>„ -> z/ T ))a,. e . 



(10) 



where k 2 has the same meaning as k\ in Eq.(8), and cr q e ,a q e are the cross- 
sections of quasi-elastic scattering of muon neutrino on (bounded) neutrons and 
anti-muon neutrino on protons in the H 2 0-target, respectively, which were taken 



to be constant and to be equal to its maximal values of a v = 1.0 x 10 38 cm 2 



q.e. 



and a qe = 0.4 x 10~ 38 cm 2 in the energy region E u <1.33 GeV, see e.g. [11 
The average cross section a q e = (a q e + a q e )/2, since according to MC simu- 
lation used in the experiment the ratio N v /N v = 1 within a few percent for 
the neutrino energy range considered ]T2|]. Thus, the total muon neutrino flux 



is underestimated. The number of detected quasi-elastic muon-neutrino events 
can be extracted from the number of /z-like events in Table 1 of ref.0. Under 
the — > v T oscillation hypothesis used the fractions of z/ e CC, u^CC and NC 
events in a detected sample of //-like events calculated for P{y^ — > u T ) = 0.3 
should be 0.7%, 94% and 5.3%, respectively, instead of 0.5%, 95.8% and 3.7% 
given in Table 1 of ref.0. Thus, N£ e - = 900 x .94 x 916.9/1166.5 = 665 events. 
We also assume that the final state signature in water Cherenkov detectors are 
practically the same for quasi-elastic u e and v T — e elastic scattering processes, so 
detection efficiencies for these processes are the same. The efficiencies for identi- 
fying quasi-elastic v e and events were 93% and 95%, respectively [0]. Finally, 
since detection efficiences for fully-contained e-like (P e >100 MeV/c) and /i-like 
(Pn > 200 MeV/c) events are practically energy independent we neglect this 
small difference in efficiencies and assume that k\ = k 2 . 
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Using the above value for N^ e -, Eqs.(7-10) and \i Vt = 5.4 x 10~ 7 fi B , which 
corresponds to the BEBC experiment upper limit on diagonal tau neutrino mag- 
netic moment [13], it is found that in the first bin of histogram in Fig. 4 (a) []J the 
S-K experiment should detect: 



AN} = m e - ■ -—4=^ ^ ~ 1420 events (11) 

This number is much greater than the number of observed events, N^ data = 
160 events, or the number of events predicted by MC simulation, N* MG = 125 
events in this energy bin. | 

The bound for the tau neutrino magnetic moment can be calculated by using 
the following relation : 



P(V„ — > Z/r 



NT ■ - - " T \ ■ ^ < (AiV, 1 ) 90 (12) 

where (AiV, 1 ) 90 = {N^ data — N^ MC ) 90 (~ 80 events) is the 90% C.L. upper limit 
for the expected number of e-like events from v T — e scattering in the first energy 
bin of Fig.4(a) (ref.JJ]), calculated taken into account the systematic uncertainty 
in the absolute normalisation of the MC events of 25% [0 by added the errors in 
quadrature to the statistical error. 

This results in a conservative bound 



fi UT < 1.3 x 10~ 7 n B (13) 
The bound is a factor 4 better than the previously published bound obtained 



by the BEBC experiment |13j and is obtained on the assumption that — > v T 
oscillations are the origin of the anomaly in atmospheric neutrino data observed 
by the S-K experiment. It can be improved by a more detailed analysis of the 
S-K data. Note that Eq.(4) and limit of Eq.(13) are consistent for the region of 
Am 2 ~ (10~ 3 - 10~ 2 ) eV 2 suggested by the S-K analysis @], E u ~ 0(1 GeV) and 
for B equal to the average magnetic field of the earth (B<1 Gauss). 

In the S-K detector, the oscillation scenario was checked by using ir° events 
IJ. For — > v T oscillations the number of neutral -current (NC) events should 
be unchanged by neutrino oscillations. For — > v T oscillations the number 
of v e CC events ( e-like events) should also be unchanged by neutrino oscilla- 
tions. Therefore, the (ir°/e) ratio of the data should agree with the same ratio 
of the Monte Carlo without oscillations for — > v T case M. It was obtained 



These numbers were read off Fig. 4(a) of ref. 
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{^/e) data /(-K Q /e)MC = 0.93 ± 0.07 (stat) ±0.19(syst) ( preliminary). The result 
is consistent with the — > v T interpretation of the S-K data, however it cannot 
exclude — > u s oscillation hypothesis completely. 

In a recent publication flip , it was pointed out that a large diagonal and/or 



transition magnetic moment of neutrino can contribute to the neutral current 
effects ( by a single 7r° production) used to distinguish the mechanisms of muon 
neutrino oscillation to tau neutrino or to a sterile neutrino. The effect discussed 
in the present Letter, as well as the effect discussed in ref. [H], definitely can 
affect the interpretation of atmospheric neutrino data. 

Finally note that helicity flipped states of neutrinos with the large magnetic 
moment would be trapped in the SN1987a core, so there is no limit on the neutrino 



magnetic moment form SN1987a |T5[. Others astrophysical constraints on /i Vr 
are also model dependent and have considerable theoretical and experimental 
uncertainties. 
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